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We consider the groundstate wave function and spectra of the L-site XXZ U q [s£(2)] invariant quantum spin chain with 
q = exp(7ri/3). This chain is related to the critical Q = 1 Potts model and exhibits c = conformal invariance. We show that 
the problem is related to Hamiltonians describing one-dimensional stochastic processes defined on a Temperley-Lieb algebra. 
The bra groundstate wave function is trivial and the ket groundstate wave function gives the probabilty distribution of the 
stationary state. The stochastic processes can be understood as interface RSOS growth models with nonlocal rates. Allowing 
defects which can hop on the interface one obtains stochastic models having the same stationary state and spectra (but not 
degeracies) as the XXZ chain. 

PACS numbers: 02.50.Ey, 11.25.Hf, 05.50.+q, 75.10.Hk 



I. INTRODUCTION 

It is known that the Q = 1 Potts model 0] is related 
to other models in statistical mechanics such as perco- 
lation Q and the O(l) loop model ||. It was only re- 
cently realized, however, that at criticality it is related to 
other topics in physics and mathematics. Consider the 
ferromagnetic quantum U q [s£(2)] invariant XXZ Hamil- 
tonian [01 



L-l 



H 



XXZ 



1 1 

4 ~ 2 



(^i+^i+ A -K+i) 



h(a*-a* +1 



(1.1) 



(1.2) 



where a x ,a y , a z are Pauli matrices, q = e 7 ™/ 3 and 

2A = q + q- 1 = ^Q = 1, h= 1 1 {q-q- 1 ) = i^ (1.3) 

The value A = 1/2 relates to the critical Q = 1 Potts 
model. The ej are the generators of a Temperley-Lieb 
(TL) algebra 



e j e j±l e j 



0, |t-j|>l. (1.4) 



The groundstate energy is zero, that is Hxxz\0) = 0|0), 
for any number of sites as can be seen by taking the 
quotient ej = 1 of the TL algebra. 

As pointed out by Read and Saleur Q the spectrum of 
the Hamiltonian (1.1) is also related, in the continuum 



limit, to the spectra of nonlinear sigma models defined 
on coset supermanifolds. Such models are relevant in 
the integer quantum Hall effect || and problems with 
quenched disorder |Q. The spectra of all of these mod- 
els are given by a generalization of the c — Virasoro 
algebras related to logarithmic conformal field theories 
(LCFTs) § (for reviews of the subject see §). What is 
common to all these theories is the Jordan cell structure 
of the representations which occur when studying their 
spectra. This is due to the fact that, at q a root of unity, 
indecomposable representations of U q [s£(2)] appear and 
the same is true for models with global supersymmetry 
(indecomposable representations are common place in su- 
peralgebras). The Jordan cell structure also appears in 
the definition of LCFTs. 

Another field where c = theories are relevant are 
stochastic models if the dynamical critical exponent z = 
1. This possibility was suggested by numerical calcu- 
lations in a 3-state model at the spinoidal point |To| |. 
In stochastic models the groundstate is again zero for 
any number of sites and, this time, the groundstate wave 
function has a direct physical meaning: it describes the 
probability distribution of the stationary state. 

In unrelated developments, study 1 1 of the XXZ chain 
has lead to various classes of alternating sign matrices 
(ASM) |l^] — a subject of deep interest in combinatorics. 
Initially the ASM numbers where found in the periodic 
chain (L odd) |h| and have subseq uently also been found 
for open |l4| and twisted boundaries pi|Jl^] . They also 
appear in the ice model with domain wall boundaries |l6| ] 
as was shown in Jl7],[l8| . This model is also related to the 
combinatorial problem of domino tiling jl9| ]. 

A relevant question is — why should number theoret- 
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ical properties of wave functions in quantum mechanics 
be significant? In this letter we answer this question by 
showing that these ASM numbers are related to directly 
measurable quantities. Specifically, we bring under one 
roof s ome of the above topics and relate the Hamilto- 
nian ( |l.l| ) to stochastic processes on TL algebras where 
the stationary state probability distribution is given by 
the groundstate wavefunction |0). A much more detailed 
version of this work will be published elsewhere |2(]] . 

II. TEMPERLEY-LIEB ALGEBRA 

The generators {ej, j — 1,2,...,L— 1} of the TL alge- 
bra can be represented graphically using monoids |2l| 



(2.1) 



3-1 j 3+1 3+2 L - 1 L 



The number of independent words Cx in the TL algebra 
with L— 1 generators is given by the Catalan numbers 

C n = — |— ( 2U ) =1,2,5,14,... 71 = 1,2,3,4... (2.2) 
n + 1 \ n ) 

We write Cj = Vjj+ivjj +1 and move to the left ideal 
by ignoring the upper half of the monoid diagrams 




12 3 4 



The vector spce generated by this ideal is equivalent to 
the state space of the O(l) loop model Q. Let us as- 
sume that L is even. Then the number of independent 
words in the left ideal is The diagrams give the 

action of the TL generators on the ideal 




e j( W J-lj+2Vij + l) = 

3-1 3 3+1 3+2 3-1 3 3+13+2 



(2.4) 
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It is convenient to encode the words in the left ideal 
by restricted solid-on-solid (RSOS or Dyck) paths \a) = 
(do, oi, . . . , cil 12) where aj is the number of half- loops 
above the midpoint between sites j and j+1 and dj+i — 
a,j = ±1 for each j. For L = 6 the RSOS paths are 



L = 6: {(0,1,2,3,2,1,0), (0,1,2,1,2,1,0), (2.6) 
(0, 1, 0, 1, 2, 1, 0), (0, 1, 2, 1, 0, 1, 0), (0, 1, 0, 1, 0, 1, 0)} 

In this basis for the TL ideal, the Cl/ 2 x matrix 
representative of H for L = 6 is 



H = -H 



xxz — 
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Since (1 — ej) = (1 — e,) 2 are projectors we see that H 
is nonpositive definite. The matrix H is representation 
independent. It is straightforward to reconstruct the cor- 
responding eigenvectors of -f/xxz using 

Vj,k = r\ = q- 1 ' 2 ] i)k - q 1/2 \ T)fc (2.8) 

j k 

where q = e™/ 3 and | |) = (1,0), | |) = (0,1) are the 
spin basis vectors. Moving to the left ideal changes the 
degeneracies in the spectra but not the groundstate wave 
function. 



III. TEMPERLEY-LIEB STOCHASTIC 
PROCESSES 

It is easy to show that for even L the Hamiltonians 
H are intensity matrices (see 0] for some properties) 
satisfying H ab > and H aa = -J2b^a H «b- The oc- 
curence of intensity matrices reveals a novel connection 
with a stochastic process with time evolution given by 
the Euclidean Schrodinger equation |23j| 



dt 



P a {t) = Y J H ab P b (t) 



(3.1) 



where P a (t) is the (unnormalized) probability of finding 
the system in the state (RSOS path) \a) at time t. In |2(]] 
we will exhibit the rules giving the transition rates H ab 
for the general case which define RSOS growth models 
(for a review of growth models see [Q). These rates 
H a b = 0,1,2 for a ^ b are simple. At each time step 
adsorption occurs with rate H ab = 1 and desorption with 
rate H ab = 1,2. These transition processes correspond to 
the addition of a single (diamond shape) tile to the RSOS 
path or the removal of a partial layer of tiles and are 
ind icat ed above and below the main diagonal respectively 
in ([j.7j). Notice that the desorption process is nonlocal. 

Since H is an intensity matrix it has a zero eigenvalue 
with a trivial bra 



(0|# = 0, 



.1) 



(3.2) 



and nontrivial ket wave functions giving the unnormal- 
ized probabilities of the unique stationary state 

ff|0) = 0, |0) =Y / P a \a), P a = lim P a (t). (3.3) 
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The wavefunction |0) describes a critical statistical sys- 
tem of weighted RSOS paths. The normalized probabil- 
ities are p a = P a /S X (L) where S X (L) = (0|0). We also 
have T <0|°1.= S ^ L ? and (°I°> T = C l/2- The conjec- 
tured Hjly] magic properties of the groundstate wave 
functions of ifxxz now assume a direct physical signifi- 
cance. 

We have calculated the wavefunctions exactly up to 
L = 18. If we now normalize the wave functions |0) to 
have smallest entry 1 



L = 6 



|0> = (1,4,5,5,11), Si(6) = 26 (3.4) 



we confirm that the normalization factors S X (L) 
(0|0) satisfy 



S x (2n) = A v (2n + 1) = 1, 3, 26, 646, . 



(3.5) 



where Ay(2n+1) is the number of vertically symmetric 
(2n+l) x (2n+l) ASMs @p| 



71-1 



A v (2n+l) = 11 (Sj + 2 



3=0 



(2j + l)!(6j + 3)! 
(4j + 2)!(4? + 3)! 



(3.6) 



We also confirm that for the highest path P a = 1 and 
for the lowest path P a — N$(L) where N$ is the num- 
ber of cyclically symmetric transpose complement plane 
partitions pljll 



N s (2n) = JJ (3j + 1) 



3=0 
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(4j)!(4j + 1)! 



(3.7) 



More generally, and to accommodate odd L, we can 
allow 2 s or fewer defects represented by unpaired (ver- 
tical) lines where s = 0, i, 1, . . . , [L/2\ is the spin. The 
number of independent words in the left ideal is then 
C(l+2s)/2- F° r example, for L = 5 



(3.8) 



1 2 3 4 5 6 



Under the action of the TL generators, the defects can 
hop and adjacent defects can be annihilated 



BjVj-xj = I C A = I ( = Vj_j + i (3.9) 

3—1 i 3+1 3 — 1 J 3 + 1 



The RSOS paths for L = 5, 2s = 1 and L = 4, 2s = 2 are 

L = 5,2a = l: {( | 0, 1,2, 1,0), (0, 1,2, 1,0, 1 1 ), 

( | 0, 1, 0, 1, 0), (0, 1, 0, 1, I ), (0, 1, I 0, 1, 0)} (3.10) 
L = 4,2s = 2: {(| 1 0, 1, 0), (| 0, 1, ]), (0, 1, 1 |), 

(0,1, 2, 1,0), (0,1, 0,1,0)} (3.11) 

The matrix representatives of H and wave functions are 



H = 



10) 
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(3.12) 



(3.13) 



The matrices are again nonpositive definite intensity ma- 
trices. For two or more defects the Hamiltonian is (Jor- 
dan cell) block triangular and the spectrum decomposes 
according to the s£(2) fusion rules which for 2 defects is 
\ <8> | = 1 © 0. For L < 17 odd and 2s = 1 we have 
verified @ that 



51(271 - 1) = N s (2n) = 1, 2, 11, 170, 7429, 



(3.14) 




It is sometimes convenient to remove the defects by 
adding 2s sites on the right, joining the defects to these 
sites by half-loops and working with extended RSOS 
paths on L + 2s sites without defects. In this scheme 
the defects are incorporated as boundary conditions. 

As for the combinatorial significance we conjec- 
ture ^0|, following [ p| , that for both even and odd L 
the integers appearing in the groundstate are given by the 
number of configurations of the fully packed loop model 
on a (2L — 1) x (\_L/2\ + 1) pyramid grid domain with 
specified boundary conditions and links determined by 
the diagrams in the left ideal. For example, for L = 6 
and v 1,2^3,4^5,6 there are 11 configurations of the type 



r~\ r~\ r~\ (3.15) 



IV. PROPERTIES OF STATIONARY STATES 

The stationary states are independent of the initial 
condition and Jordan cell stucture but the time evolution 
does depend on which Jordan blocks are selected out by 
the initial t — probabilities. A stochastic process can 
also be defined by the action of H on the whole TL al- 
gebra (Cx states) giving a different Jordan cell structure 
and (non-unique) stationary states. 

Many properties of the stationary states can be stud- 
ied by finite-size scaling. We give preliminary results on 
the average perimeter (1(a)), average area (N(a)) and 
average number of clusters (C(a)) where we define 

2[i/2]-l 2TL/2] 

1(a) = ?\"W-«J-i\, N(a)=J2 [ aj /2\ (4.1) 
i=i j=i 

L 

C(a)='£S(a j ,0)-2s, (...)= (4.2) 

j=l a 

Here N(a) is the number of tiles and we have ignored the 
fixed horizontal contribution 2\L/2~\ to the perimeter. 
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By extrapolating results for L < 18 we find numerically 
the following scaling behaviour for large even L 

(£(a)) ~ .249(1) L, (N(a)) ~ .065(1) L 1+v , v = .50(3) (4.3) 



(C(a)} ~ 1.17(1) L 1 , as = .667(3) 



(4.4) [1 



The same exponents but different amplitudes are ob- 
tained for (£(a)) and (N(a)) with odd L. The finite- 
size behaviour of other properties describing the critical 
weighted RSOS paths model and its connection to con- 
formal invariance will be discussed in 20 1. 



V. CONFORMALLY INVARIANT SPECTRA 

We assert that the spectra of the intensity matrices H 
are described by a conformal field theory |j] with central 
charge c = 0, conformal weights 



s(2s-l) 
3 



0,0 i 1 



(5.1) 



and Virasoro characters with q the modular parameter 

oo 

X s(q) = <z As (l - q 2s+1 ) J] (l - g")" 1 (5.2) 



Finite-size corrections 25 to the energy levels E n , 
H\n) = —E n \n), for large L with n = 0, 1, 2, . . . are 



LE n 



= A s + fc„ + o(l), fc„eN 



(5.3) 



where fc„ labels descendents and v — 3V3/2 [^6| is 
the sound velocity. We have calculated numerically the 
finite-size spectra up to size L = 16 and obtain the fol- 
lowing estimates using Van den Broeck- Schwartz [^7j ap- 
proximants. 
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TABLE I. Table of energy level estimates for 2s = 0, 1, 2 
defects giving the characters Xs(o) f° r s = 0, ^, 1. 
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